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SUMMARY 


A theoretical investigation is presented of the time-temperature- 
dependent buckling of a pin-jointed constant-section column, whose initial 
curvature is defined by a half-sine wave when the material is linearly 
viscoelastic and is heated uniformly along the column at a prescribed time 
rate. It was .found that the deviations from straightness increase with 
time and become indefinitely large when heating reduces the Young’s modu- 
lus of the material to the value at which the applied load is the Euler 
load of the column. When the column is heated very rapidly this critical 
time represents the limit of usefulness of the column. When heating takes 
place less rapidly the deflections of the column cause bending stresses 
exceeding the yield stress of the material at a time considerably smaller 
than the critical time. The equations presented permit the calculation 
of this reduction in the useful lifetime of the column. 


INTRODUCTION 


Analytical investigations and experiments have shown that a visco- 
•elastic column at constant temperature subjected to a constant end load 
less than the Euler load will buckle if the load is maintained for a 
sufficiently long period of time. If the temperature is increased the 
time required for buckling decreases (refs. 1 to 3)* This type of response 
of a structure to a sustained constant load is an example of the effect 
of creep, which is one of the factors responsible for the inelastic 
behavior of a col umn - 

Creep buckling has been investigated recently in some detail by 
Kempner (refs, 4 and 5)- In reference 4 the column was assumed to possess 
ideal linear viscosity of the Newtonian type. It was found that initial 
slight deviations of the center line of: the column from the straight line 
increased continuously with time and became indefinitely large if the load 
and the temperature were maintained constant for an indefinitely long time. 
The mathematical concept of indefinitely large displacements is equivalent 
to the practical concept of buckling. However,' a column becomes useless 
for practical purposes at an earlier time, namely when it becomes curved 
to such a degree that it cannot fulfill any more its structural purposes 
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or when the "bending stresses caused "by large deviations from straightness 
cannot he supported by the material of the column. The situation is 
different when the material exhibits nonlinear viscosity. In that case, 
investigated by Kempner in reference 5, indefinitely large deformations 
are reached according to the equations derived in a finite rather than 
in an indefinitely large value of time. Nonlinear ly viscoelastic columns 
buckle, therefore, with a snap action. 

In all the work cited the column was assumed to be subjected to a 
constant load wlpLle the temperature was maintained at a constant level. 

In reality the, temperature of the structure of a supersonic plane or 
guided missile is Increased from the temperature prevailing at the air- 
port or in the hangar to some elevated temperature corresponding to 
thermal equilibrium at full supersonic speed while the plane or missile 
is under the action of the loads existing in flight. It is of interest, 
therefore, to examine the problem of the buckling of a column subjected 
to a given constant load while its temperature is raised from ambient 
temperature on the ground to the thermal equilibrium temperature in super- 
sonic flight. In the present report such an investigation Is carried out 
with the assumption that heating takes place at a prescribed rate. 

The inelastic behavior of a material was analyzed by the theory of 
mechanical models whereby the constituent phases of the material were 
replaced by mechanical models (ref. 6). The mechanical models are combi- 
nations of two model elements representing the two basic types of defor- 
mation: 

(1) A perfectly elastic spring, which obeys Hooke's law, for elastic 
deformation. 

(2) A dashpot, consisting of a perforated piston moving in a cylinder 
containing a viscous liquid, for viscous deformation. 

The spring element is a model of a linearly elastic body while, if the liquid 
in the cylinder obeys Newton's law of viscosity, the aashpot is an example 
of a linearly viscous body. These basic elements -may be coupled in series 
or parallel. When combined in series, a linear Maxwell model results and 
represents the model considered in this report to explain the Inelastic 
behavior of the material. 

The total strain o.f a Maxwell model under a uniaxial force consists 
of an elastic component (which contains Young's modulus) and a creep 
component (which contains a viscosity coefficient). In this report, 

Young's modulus and the viscosity coefficient are considered to be func- 
tions of temperature, while the temperature is a function of time. Inertia 
effects are neglected, since it is shown in reference 7 that there Is good 
agreement between results of the dynamic and static approach for linearly 
viscoelastic columns. 
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Based on the assumption that during bending any transverse section, 
originally plane, remains plane and normal to the longitudinal fibers of 
the beam, the veil-known fundamental relation between radius of curvature 
and strain is obtained. Ass umi ng small deflections, the latter equation 
when differentiated with respect to time results in a third-order, linear 
partial differential equation. The assumption of a product solution in 
the independent variables reduces this equation to a first-order equation 
whose variables are separable. 

The deflection of the column is expressed in terms of a simple, 
definite time integral containing Young's modulus and the viscosity coef- 
ficient. A solution is presented for a general linear variation of Young's 
modulus and viscosity coefficient with temperature and of temperature 
with time.' 

The author wishes to express his gratitude to Professor Joseph 
Kempner for his advice and to Professor IT. J. Hoff, the supervisor of 
the project, for suggesting the problem. 

This work was conducted at the Polytechnic Institute of Brooklyn 
under the sponsorship and with the financial assistance of the National 
Advisory Committee for Aeronautics. 


SYMBOLS 


A 

a l> b l 


®3 , b 3 


C 

E 

E, 

F 


cross-sectional area of column 

constants; coefficients in assumed linear Young's modulus - 
temperature relation 

constants; coefficients in assumed linear viscosity-coefficient - 
temperature relation 

constants; coefficients in assumed linear temperature-time 
relation 

constant, EoP/Pe q 

Young's modulus 
activation energy 
defined by equation ( 22 ) 
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maximum deflection of unloaded column 
I moment of inertia with respect to neutral axis 

•K constant defined by equation (19) 

L length of column 

M moment of external forces 

P constant, compressive end load 

P E (t) temperature (time) dependent Euler load, ir^E(t )l/L^ 

R universal gas constant 

T temperature 

T_ temperature, deg abs 

. t time 

w additional deflection from initially curved, unloaded column 

v i initial deviation from straightness 

v 0 elastic deflection resulting from instantaneous application of 

end load P 

x axial coordinate of column 

z normal distance from neutral surface 

a constant, a 2 /h 2 

3 constant, (aq - C)/b^ 

e total axial strain 

e c creep component of axial strain 

e e elastic component of axial strain 

X viscosity coefficient 

p radius of curvature 

a stress 
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4> defined by equation ( 13 ) 
t£- function of temperature 
Subscript: 

o conditions at zero time 


GENERAL ANALYSIS 

The fundamental "bending relation between the total strain 6 and 
the radius of curvature p, which is based on the assumption that any 
transverse section, originally plane, remains plane and normal to the 
longitudinal fibers of the beam after bending, is 

l/p = e/z (l) 

where z is the normal distance from the neutral surface. 


The inelastic behavior of the column is expressed in terms of a 
linear' Maxwell unit, that is, a series combination of a spring element 
obeying Hooke’s law for linear elasticity and a dashpot whose fluid obeys 
Newton’s law of viscosity for linear viscosity. The total strain of the 
Maxwell unit under a uniaxial force consists of the elastic component 


and the creep component e c (ref. 6), or 

= (cr/E) + J ( a/X ) dt 


e = 6 e + e c 


( 2 ) 


where Young 1 s modulus E and the viscosity coefficient X are both 
functions of temperature, while the temperature is a function of time t. 


The substitution of equation (2) into equation (l) and multiplication 


t>y 


yield 


Therefore 


z 2 /p = (az/E) + z f (a/X) dt 


z 2 d(l/p)/dt = (ctz/e) - (Ecrz/E 2 ) + (az/X) 


where d/dt = (*). Integration of each term over the cross-sectional 
area A results in 


Afi/p^dA 


(l/E) J h 


az dA 



crz dA + 



az dA 


( 3 ) 
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As the moment M of the external forces is expressed by 


*- 1 


a z dA 


it follows that 


U.fo 


a z dA 


Hence equation ( 3 ) reduces to 

El d( l/ p)/dt = M + [(e/X) - (E/E)]M (4) 

where I is the moment of inertia with respect to the neutral axis . 

For small deflections 

l/p = -d^w/ dx 2 (5) 


where w is the additional deflection from the initially curved, unloaded 
column and x is the axial coordinate (fig. l) . The initial displace- 
ment w.^ of the pin- jointed, unloaded bar is assumed to be 


W£ = sin (rtx/L) 


( 6 ) 


while the elastic deflection w Q resulting from the instantaneous appli- 
cation of the constant, compressive end load P is obtained from the 
solution of 

E q I d^/dx 2 = -Mo = -P(w Q + w ± ) (7) 


The solution of equation (7) is 


P/P E 

w o = F i 7 — ~ — \ sin (Wl) 

1 - ( p / p O 


(8) 


where the subscript o refers to conditions at t = 0, and Pw = 

u 

is the corresponding Euler load. When t > 0, the creep phenomenon must 
be considered) then the external moment M at any section is 


M = P(w^ + w) 


(9) 


so that- 


M = Pw 


( 10 ) 
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The substitution of equations (5), (9), and (10) into equation (4) 
and the rearrangement of terms result in the following t hir d-order, 
linear partial differential equation: 


c^w 


EI + Pw + P[(EA) - (E/E)]w = -P[(EA) - (E/E)] Wi 


Equation (ll) must satisfy the conditions: 


( 11 ) 


t = 0, w = w 0 
x = 0, L, w = 0 


As a solution, assume 

w = F(t) sin (jtx/l) (12) 

which, when introduced into equation (ll), together with equation (6), 
gives 

-EIFU/l) 2 + PF + p[(eA) - (e/e)]e = -p[(eA) - (e/e)^ 


Let 


P E = Pe^) = «%(t)l/L 2 


0 = <S(t) = [(eA) - (e/e)] /(E - C) 


(13) 


where C 


E 0 P/1 


E, 


is a constant; therefore the last equation becomes 


F - C$F = C$Fj_ 


Hence 


whose solution is 


dF/(F + F ± ) = C<t dt 



- E ± 


F = (F 0 + Fj.)e 
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From equations (8) and (12), when t = 0 


therefore 


w = w, 


= F. 


P A 


1 1 - ( p A*o) 
= F sin ( itx/L) 


sin (joc/l) 


F 0 = Fi 


P/P. 


= FiC/(Eo - C) 


and 


F 


?A 


[V< E ° - 




J 


cut} 


This expression is further simplified, by the evaluation of the 
integral 

r\ dt . r tvy)_- at 

J 0 J 0 (E - C) 

nt r>E 

-Jo I(F 1 F) at -J Eo [dE/E(E - c,] 

= jC at - (1/c) T{ [1/(E - c)] - (1/E) ) “ 


— 

1 0 X.(E - c) 


dt - ( l/C) log e [(E - C)E 0 /(E 0 - C)e} 
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Therefore 


CjT ® dt CJ^ __S_^dt-log e [(E-C)E 0 /(E 0 -C)E] 


= e 


E 


= [e(E q - C)/E 0 (E - c)] e 
Equation (l4) becomes 


! 0 X(E-C) 


dt 


F = E, 


of 

J 0 


E 


X(E-C) 


-dt 


- 1 


1 - (C/E) 


But C/e = E 0 P/Pg E = P/Pgj hence 


F = 




E 


>0 x(e-c) 


dt 


1 - (P/Pfi) 


- 1 


(15) 


and from equation (12) 


v = Fi 


C / 

^ r\ 


E 


0 X(S-O) 


1 - ( P / P E> 


dt 


- 1 


sin (jdc/l) 


(16) 


Equation (l6) expresses the additional deflection w in terms of a 
simple, definite time integral containing Young's modulus E(t) and the 
viscosity coefficient X(t). It will be noted that the deflected column 
retains its sine-curve characteristic and that the additional deflection 
tends to infinity In a finite time when the temperature (time) dependent 
Euler load P-g(t) approaches the applied, constant end load P. 

Although the equation indicates that the infinite deflection is inde- 
pendent of the viscosity coefficient, the time -dependent viscosity coef- 
ficient does, for practical purposes, influence the buckling failure if 
buckling failure is defined by some finite ratio of additional deflection 
to initial deflection. 
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APPLICATION OF GENERAL DISPLACEMENT EQUATION 


Equation ( 1 6 ) is solved for the following assumed general linear 
relations between Young's modulus and temperature, viscosity coefficient 
and temperature, and temperature and time: 

E = a-j_ + bjT 

X = a 2 + b 2 T 

T = a^ + b^t 

r 5 dt = C f [dt/X(E - C)] + f (dt/X) 

J 0 X(E - C) J 0 

= (cAlV^) / [<3T/(T + a)(T + p)] + 

J T 0 

(i/b 2 b 3 [ar/(T + a)] 


where a - P = ( a i ” < -')/' 3 l* 


^ [ffC/(T + a)(T + p)] = -[l/(a - p)] f [®/(T + a)] + [l/(a - p)l J [d3?/(T + p)] 
J T 0 J T 0 T o 

= - [l/(a - p)j loge (T + a)j + j[l/(a " p )] lo 6 e + 


T 

4 T 0 


. [l/(a - fO] lo fe ; t ± si 

(T 0 + P)(T + a) 


= log- 


f (T 0 + a)(T + 
j_(T + a)(T 0 + 


P) 

P)| 


|a-P 
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and 


r i> /( t + a >i = logg [(T + a)/(T 0 + a)] 


therefore 


r 


0 uFTst dt ' ( c/b iV 3 ) ^ 


(To + a)(T + p) 


1 _ 
rra-p 


(T + a)(T 0 + p) 
(l/bg^) logg JjT + a)/(T 0 + a)] 


Let 


= lo Se 1 


(T 0 + a)(T + p) 


(T + a)(T c + p) 

• — 



T + a \ bgb 3 1 


/ T + a \ 

Uo + 


f H 



C 1 1 

(T 0 + a)(T + p) 

T^Ca-p) t + ,^3 | 

|_(T + a)(T 0 + p) 

Uo + a/ ■ 

•w 

J 
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Then 


c f*. 

J 0 X(E-C) 


C logelr 
e 


. e lo 6e 


= * 


C 


C 


2 


(Tp + a)(T + (3) 
(T + <x)(T 0 + p) 





. . c g . c 2 

/ T 0 + q \ b l b 2 b 3 (a " p ) V3 / T + p \ b l b 2 b 3^ a -P) 

\ T +cj \T 0 + 3/ 


(17) 


Reference 8 presents an empirical quadratic formula for Young’s 
modulus variation with temperature for 75S-T6 aluminum. The quadratic 
formula is assumed to he approximated by the following linear relation: 

E = (l0.5 X 10 6 ) - (7.5 X 10 3 t) (18) 

where T . is in. degrees Fahrenheit. 


In general, the viscosity-coefficient variation with temperature 
(ref. 9) Is expected to be of the form 


Jk 

\ = Ke^a 


( 19 ) 


where 

E q activation energy 

R universal gas constant 

T a . temperature, deg abs 


K 


constant 
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The approximation of a linear relation to an exponential curve is 
valid only for a narrow range of the temperatures. Based on the experi- 
mental curves of figure. 3 in reference 2 and the use of equations (2) 
and (19), the following relation “between viscosity coefficient and 
temperature is assumed: 


X = 


(6 x 10 13 ) - ( 1.5 X 10 11 t) 


( 20 ) 


It is assumed that at zero time the temperature of the column is 
75° F and that its temperature rises to 400° F within 60 seconds. There- 
fore the temperature-time equation is 

T = 75 + 5.42t (21) 

where t is the time in seconds. 


A comparison of equations (18), (20), and (21) with the general 
linear equations gives 


a-L = 10.5 X 10 c 
ag = 6 x 10 13 


a 3 


b ± = -7.5 x 10 J 


b 2 = -1.5 X 10 

b 3 = 5.42 


11 


If P/P E =0.8 


c = e q p/p Eo = 9.9375 X 10 6 X 0.8 = 7.95 X 10 6 
a = B^jb^ = -400 
P = ( a i - c/h-L = -340 
C 2 C 


b-jb^b^a, - f3) b2b^ 


= -O.OOOI63 


7 = -0.000173 

b l b 2 l 3 3 ^ a * 
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Therefore equation ( 17 ) becomes 


E 


I 0 X(B-C) 


dt 


(Tp + p)°- 0001 73( T t al O.O°Ql63 
(T 0 + o) 0.000J6 3(i + p)0. 000173 


T - 4oo| °< °0° 163 
T . 3 4 o| 0.000173 


since 


<T 0 + S ) 0 ' 000173 

75 . 3 to |o.oooi73 

(t 0 + a ) 0 ' 000163 

73 - tool 0 ' 000163 


Hence equation ( l6) is 

| T . tool 0 ' 000163 


W 

P. 


1 [l - (C/E)]| T - 34o|°' 0001T3 


.> sin (kx/l) 


( 22 ) 


This relation is shown in figure 2 for the midpoint deflections of 
the column. Also shown in figure 2 is a curve for a column subjected to 
the following conditions: 



0.9 


and 


T = 75 + 0.5U2t 

Of the two curves the steeper one corresponds to a rapid heating 
with a critical time of 48 seconds. The rate of heating corresponding 
to the flatter curve raises the temperature of the column in such a manner 
that the Euler load is reached in 244 seconds . Consequently, even the 
slow heating is rapid enough from the viewpoint of practical applications. 

The diagram clearly indicates that indefinitely large deflections 
develop at the time when the increase in temperature lowers the Young's 
modulus of the material to such an extent that the applied load becomes 
the Euler load of the column. At the more rapid rate of heating presented 
the critical time can be considered in good approximation as the limit of 
practical usefulness of the column. On the other hand with the slow rate 
of heating the column becomes useless at a value of the time considerably 
smaller than the critical time. The cumulative effect of creep results 
in deflections of such magnitude that the linear stress -strain relation- 
ship loses its validity; the material starts to yield rapidly and the 
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column truckles for all practical purposes before the theoretical critical 
time, based on idealized conditions, is reached. If a maximum deflection 
of the column amounting to 50 times the initial maximum, deflection is 
considered inadmissible at the slow rate of heating, the column would 
become useless after 200 seconds rather than at the critical time of 
2kb seconds. 


DISCUSSION 


The general displacement equation (eq. (l6)) indicates that the 
deflected column retains its sine-curve characteristic and that the addi- 
tional deflection tends to infinity in a finite time when the temperature- 
dependent Euler load approaches the applied, constant end load. 

Although the equation states that the infinite deflection is inde- 
pendent of the viscosity coefficient, the time- dependent viscosity coef- 
ficient does, for practical purposes, influence the buckling failure if 
buckling failure is defined by some finite ratio of additional deflection 
to initial deflection. This can be seen from figure 2 which gives the 
midpoint nondimensional additional deflections as a function of time for 
the assumed linear relations. At the more rapid rate of heating the time 
required for the temperature-dependent Euler load, to attain the value of 
the applied constant end load is U8.8 seconds. With the slow rate of 
heating this critical time is 2 W- seconds. If 50 times the initial maxi- 
mum deflection is considered as the limit of structural usefulness of the 
column the lifetime of the column is 200 seconds and not the theoretical 
critical time of 2kh seconds. 


Polytechnic Institute of Brooklyn, 

Brooklyn, N. Y., July 23, 1952. 
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